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Fragmentation reactions with intermediate-energy heavy-ion beams exhibit a wide range of reac- 
tion mechanisms, ranging from direct reactions to statistical processes. We examine this transition 
by measuring the relative population of excited states in several sd-shell nuclei produced by frag- 
mentation with the number of removed nucleons ranging from two to sixteen. The two-nucleon 
removal is consistent with a non-dissipative process whereas the removal of more than five nucleons 
appears to be mainly statistical. 



PACS numbers: 25.60.-t,25.70.Mn,27.30.+t 

Beams of exotic nuclei provide a key tool to study nu- 
clear structure far from stability. At intermediate en- 
ergies (~100 MeV/nucleon), exotic beams can be pro- 
duced via in-flight fragmentation, where the nucleus of 
interest is obtained by removing nucleons from an inci- 
dent primary beam. During this reaction, the population 
of excited states depends on the number of nucleons re- 
moved. One-nucleon removal reactions from projectiles 
at intermediate energies have been shown to be mainly 
direct processes with the resulting cross sections to 
different excited states sensitive to the quantum num- 
bers of the removed nucleon. Two-nucleon knockout has 
recently been shown to be direct in certain energetically 
favorable cases 0, 01 ■ Fragmentation of heavy nuclei, 
on the other hand, can be described as a two-step pro- 
cess commonly called abrasion-ablation During the 
abrasion stage the nucleons in the overlap volume of pro- 
jectile and target are scraped off as the ions pass each 
other. In the ablation step the prefragment dissipates its 
excitation energy gained during the abrasion by particle 
emission |5j and the different states of the residues are 
populated statistically with weights determined by the 
excitation energy of the prefragment. Few- nucleon re- 
moval reactions lie between these two limiting scenarios 
of purely direct and statistical production. Discrepan- 
cies between observed and predicted excited state pop- 
ulations fell indicate the presence of non-statistical 
processes in few-nucleon removal reactions. In this ar- 
ticle, we present new data on the relative population of 
bound excited states produced in fragmentation reactions 
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and show that several-nucleon removal reactions at inter- 
mediate energies can be understood in terms of a smooth 
transition from direct to statistical processes. 

In the experiment, we measured 7-rays in coincidence 
with fragmentation residues to determine the population 
of each of their bound states. Fragmentation residues 
were obtained from ■^'^Ar, ^^Si, ^^Al and ^^Mg secondary 
beams at ~105 MeV/nucleon, which in turn were made 
by fragmentation of a 150 MeV/nucleon '^^Ar primary 
beam from the Coupled Cyclotron Facility at the Na- 
tional Superconducting Cyclotron Laboratory (NSCL) 
at Michigan State University. The ^Be production tar- 
get was located at the mid-acceptance target position of 
the A1900 fragment separator 9]. The exotic secondary 
beams impinged on a 376 mg/cm^ ^Be secondary tar- 
get located at the target position of the high-resolution, 
large-acceptance S800 spectrometer The magnetic 

rigidity of the spectrometer was set to transmit the 
residues of interest, and prevent direct beam from en- 
tering the focal plane. Each incident beam particle was 
characterized event-by-event in a time-of-flight measure- 
ment. Reaction residues were identified in the S800 focal 
plane through measurement of the energy loss in the ion- 
ization chamber of the spectrometer versus the time-of- 
flight measured between the object point and the focal 
plane. Gamma-rays from residues produced in bound 
excited states were detected in coincidence to determine 
the population of excited states. The secondary target 
was surrounded by the SeGA array |llj | composed of sev- 
enteen 32-fold segmented, high-purity Germanium detec- 
tors positioned at 20.8 cm from the center of the target. 
Seven detectors were placed at forward angles in a 37° 
ring with respect to the beam axis and ten detectors were 
located in a 90° ring around the target position. The pho- 
topeak efficiency of ScGA, in this configuration, is 2.5% 
for a 1.3 MeV photon emitted from a moving source with 
a velocity of v = 0.35c. The response of the detector was 
simulated with the GEANT code |3 a-nd the simulated 
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FIG. 2: Level-scheme of ^^Na. Experimental data (left) are 
compared to shell- model calculations (right). 



FIG. 1: 7-ray spectrum in coincidence with Na residues 
produced from ^''Mg. The data points are shown with statis- 
tical error bars. The individual components of the simulated 
spectrum corresponding to individual 7-rays, the exponential 
background, and their sum (thick solid line) are presented. 

detection efficiency validated with calibrated-source mea- 
surements. The population of excited states was deter- 
mined from the 7-ray intensities relative to the number 
of reaction residues. The population to the ground-state 
was evaluated by subtraction. More than 90 % of the 
momentum distribution was transmited for the studied 
two-nucleon removal reactions. For the other reaction 
channels, only part of the momentum distribution was 
within the acceptance of the S800, therefore no absolute 
cross sections could be measured. In the following, we 
assume that the momentum acceptance does not affect 
our conclusions. 

The fragmentation residues analyzed here were chosen 
(i) to have a separation energy of less than 5 MeV to 
exclude potentially undetected high-energy feeding tran- 
sitions and to thus ensure reliable relative population of 
individual excited states, and (ii) to be experimentally 
well-known and reproduced by shell- model calculations. 
With these criteria, we focused on the production of the 
sd-shell nuclei ^^Na, ^^Ne and ^"^Si, produced from dif- 
ferent incoming beams. 

^^Na (Sp= 2.43 MeV) was produced from fragmenta- 
tion of ^"^Mg, ^^Al and ^''Ar. Its 7-ray spectrum in coin- 
cidence with incoming ^"^Mg is shown in Fig ^ -^^Na 
presents three bound excited states and an unbound 
state that decays via electromagnetic transitions (T^ 3> 
Tp 01 )• The level scheme of ^-"^Na is compared to shell- 
model calculations performed with the code Oxbash 
and the USD interaction il^ i n Fig.|21 The 7/2+ excited 
state decays with 93(2)% [ig to the 332 keV low-lying 
level via a 1384 keV transition, whereas the 2424 keV 
state decays directly to the ground-state. The proton- 
unbound 2829 keV level is known to decay via 7 emis- 
sion with a 63(5)% branching ratio to the 7/2+ excited 
state and with 37(5)% to the 332 keV level. The rela- 
tive populations of these states in fragmentation of ^"^Mg 
and '^'^Ar are presented in Fig. 13 A first comparison of 
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FIG. 3: Relative population of excited states in ^^Na from 
^*Mg and ^*Ar. For each reaction, the predictions PI"" for 
a non-dissipative process and the fit with the population 
Pfc=Pr+Pfc are shown. 



the experimental data indicates that the fragmentation 
from '^''Ar favors the high-spin states 7/2+ and 9/2+ rel- 
ative to the fragmentation of ^''Mg. The theoretical lines 
are discussed later. -'^^Ne was produced from ^"'Mg, ^^Al, 
^^Si and ■^^Ar. It has four bound excited states below its 
proton separation energy of Sp=3922 keV. A 2+ excited 
state at 1887 keV and a (0+,2+,4+) triplet at 3500 keV. 
The members of the triplet decay to the first 2+ state. 
The relative population of states in ^*Ne from the frag- 
mentation of '^^ Ar is presented in Fig. 0] The production 
of ^''Si from ^^Si was also studied as an interesting case 
of two-neutron removal. Two transitions were observed 
in ^"^Si (Sp= 3.3(4) keV) corresponding to two excited 
states previously reported in [13 • They were assigned 
as the first and second 2+ states of ^^Si at 1860(8) keV 
and 3410(14) keV respectively, based on good agreement 
with both shell-model calculations and the level scheme 
of the mirror nucleus '^^Ne. 
In the following, we present a simple approach to qual- 
itatively reproduce our observations. The main idea is 
to consider the coexistence of both, statistical feeding of 
different states and a non-dissipative process which de- 
pends on the structure of the projectile and the residue. 
In a first step, we consider a statistical population. As- 
suming a two-step fragmentation process with a ther- 
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FIG. 4: Relative population of excited states of Ne from 
^^Ar. The predictions from non-dissipative population Pfe™ 
(dotted-dashed line) , statistical population with a temper- 
ature T=l MeV (dotted line) and T=3 MeV (dashed line), 
and best fit from Eq.|Sl (solid line). 



malized prefragment and neglecting the effect of particle 
emission on the relative feeding of excited states in the 
residue, the statistical population of a residue's bound 
excited state fc, with a spin and an excitation energy 
E^, is expected to be proportional to the Boltzmann fac- 
tor weighted by the spin degeneracy 

= N{2ju + l)e-^^*/T^ (1) 

where N = l/(X]fePfc) f^i' normalization. In this 
parametrization, T is not the usual temperature of the 
prefragment since P| is an estimation function for the 
population of the excited states of the residue. T rep- 
resents an average quantity including the effects of the 
temperature of the prefragment and the different decay 
channels of the prefragment to the residue. P** favors 
the population of high-spin and low-energy states. These 
statistical predictions are plotted for ^*Ne in Fig. 21 for 
r=l MeV and r=3 MeV. For the lower T, the popula- 
tion decays exponentially with the excitation energy E*. 
At r=3 MeV ^ Sp=3.9 MeV, the population is driven 
by both the spin degeneracy of the excited states and the 
exponential decay with E*. Clearly, no value of T can re- 
produce the observations within error bars. 

In a second step, we consider the population of ex- 
cited states to be only due to non-dissipative processes, 
such as direct multi-nucleon removal or sequential strip- 
ping reactions. To estimate the relative populations of 
such mechanisms in a simple picture, we assume (i) that 
the relative production in these non-dissipative channels 
depends only on the projectile and residue shell-model 
wave functions and the single-particle removal cross sec- 
tions, and (ii) that the cross sections for the removal of 
an s-wave and d-wave nucleon from the sd-shell are iden- 
tical, as it is qualitatively justified by the single particle 
cross sections reported in We estimated an effective 
multi-nucleon removal strength between the projectile 
and the residue via shell-model wave functions. In a shell- 
model description, the wave- function |$o) of the incom- 
ing particle with A nucleons is decomposed into its un- 
derlying independent-particle shell-model configurations 
0j(A)) with weight coefficients = |((?!)i(A)|$o)P- In 
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FIG. 5: Relative production of excited states of Ne from 
^®Mg Experimental data are compared to theoretical pre- 
dictions, namely the population calculated from Eq. |31(Pfc™), 
a complete stripping calculation (knockout) and a statistical 
population (PJ) with T=l MeV and T=3 MeV. 

the same way, the wave-function \^k) of the residue, 
produced in an excited state k from a AA-nucleon re- 
moval, is decomposed with the weight coefficients (3^ = 
{(pj (A—AA) 1 P . The effective transition strengths Ck 
between the incident particle and the residue in a state 
k, resulting from the removal of neutrons and protons, is 
taken as 

i j 

where Cy is the number of ways one can remove AA 
nucleons from \<j>i{A)) to obtain — AA)). is 

evaluated as 

c,^=l[C{nt{A),n'^{A-AA)), (3) 



where C{nf, n^) is the number of combinations of nu- 
cleons choosen from a set of nucleons, and n^{A) is 
the number of nucleons of the configuration \(j)i(A)) in 
the subshell 9 = {£,j,T), with £, j and r the angular 
momentum, total angular momentum and isospin of the 
considered subshell, respectively, is non-zero only 
if the number of nucleons n^{A) in each occupied or- 
bital 6 of \(f>i{A)) is greater or equal to the corresponding 
number of nucleons n^j{A — A A) in the same orbital of 
|(/)j(A — AA)). For simplicity, angular momentum cou- 
plings are neglected. The shell-model prediction P|™ for 
the population of each bound state k is then given by 

pr = (4) 

The calculation of P|"' for a given system is parame- 
ter free. All calculations were performed within the sd 
shell-model space with the shell-model code Oxbash 
and the USD interaction The difference between 

P^.™ and a complete calculation is shown in the case of 
^^Ne produced by the two-proton knockout from ^*Mg 
in Fig. [S] The stripping results, considering exact two- 
nucleon amplitudes and reaction cross sections, are taken 
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from ref . 's'l . Eq. 0] reproduces the trend of the complete 
calculation which is in good agreement with the data. 
This two-nucleon removal has been shown to be a direct 
process 0. Thus a standard direct reaction mechanism 
is expected to predict the experimental feeding correctly. 
In this case, the statistical population (from Eq. ^ can- 
not reproduce the experimental behavior of the relative 
population within the experimental uncertainties, espe- 
cially the strong population of both the ground state and 
the 4+ state at E^=3500 keV. The dependence of P"™ 
on the structure of the projectile for a given residue is 
clearly apparent in the case of ^^Na produced from ^^Mg 
and ^''Ar (Fig. I^J. For a ^*Ar projectile, the population 
function P'*™ is fiat, whereas it shows large differences 
from state to state for incoming ^''Mg. 
Neither the statistical assumption nor the derived equa- 
tion Eq. ^ can reproduce the entire data set correctly. 
This is not surprising since the number of removed nu- 
cleons ranges from AA=2, known to have a dominant 
direct part in the production mechanism, to AA=16 ex- 
pected to be well described by a two-step fragmentation 
statistical calculation. 

For this reason, we considered the coexistence of the 
two phenomena. The relative populations of each of the 
nine aforementioned systems were fitted with a function 
Pfc combining the statistical feeding P^ with the non- 
dissipative population P^™ in the following way 

P,=aPr + (l-«)P^(T), (5) 

where the weight of the single-particle component a 
and the temperature-like parameter T of the statistical 
part were the only parameters. Both P'' and P^™ are 
normalized to unity so that J2k-^k = 1- In adding the 
two probabilities, we neglected interferences between 
the two processes. For the two-nucleon removal, we 
considered the complete stripping reaction calculation 
for the direct single-particle function p^™. The fitting 
was performed by minimizing a standard function 
with the following boundary conditions: < a < 1 and 
< T. The results obtained are summarized in Table 
For each system, we report the values for the fit with 
Eq. [Sl(Pfc), for the non-dissipative estimation P^™, and 
for the fit with the statistical part P|. The different 
values are not normalized to the number of degrees of 
freedom (depending on the number of parameters in the 
considered probability law) in order to offer an unbiased 
evaluation of the agreement of the data set and the 
theoretical assumption. The best fit with P^ for each 
system reproduces the data well, as is illustrated by the 
good values. For the residue ^^Ne, the experimental 
uncertainties corresponding to the ground state and the 
first excited state are larger relative to other residues 
due to the decay scheme of ^^Ne. Therefore, the fit 
is less sensitive, resulting in large error bars for the 
estimation of a. The fit for the two-proton removal 
9Be(28Mg,26Ne)X gives a = 0.97(11), in agreement 
with Q. For cases with the residue corresponding to 
AA=6-16, the fit is consistent with a = 0, showing 



TABLE I: Comparison of data from this work and 2] with 
different assumptions (see text). The different reactions are 
classified with increasing number of removed nucleons AA. 
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FIG. 6: Correlation between the weight a of direct process 
(see text) and the temperature-like parameter T of the sta- 
tistical part. 



that for AA > 5 statistical population can account for 
the observations. Even if not strictly monotonic, the 
decrease of a is clearly correlated with the number of re- 
moved nucleons AA. The weight a of the non-dissipative 
removal and the temperature-like parameter T of the 
statistical process show a strong correlation as illustrated 
in Fig. From the simple ingredients of the model, 
the observed correlation can be intuitively understood. 
At T=0 the population is expected to be non-statistical 
since P''=0 except for the ground-state, whereas for 
increasing T, the stastistical feeding should be more 
dominant. This correlation validates the underlying 
intuitive assumptions of our approach. 

We have fragmented sc?-shell nuclei at intermediate 
energies and measured the relative feeding of excited 
states of the residues via 7-ray spectroscopy. Nine 
multi-nucleon removal reactions from AA=2 to 16 have 
been studied. The two-nucleon removal is confirmed to 
be consistent with a non-dissipative, direct mechanism 
whereas the removal of more than five nucleons is 
consistent with a purely statistical process. The results 
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confirm that both a non-dissipativc process and a sta- 
tistical mechanism participate in multi-nucleon removal 
reactions when 2 < AA < 6 are being removed. 
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